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state atom placed within an arbitrary dispersing and absorbing magnetodielectric multilayer system 
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I. INTRODUCTION 

It is well known that an atom near a neutral macro- 
scopic body is subject to a force, even if the atom and 
the body are in the (unpolarized) ground states. The 
existence of the force commonly called van der Waals 
(vdW) force has been experimentally well established. 
In particular, vdW forces on atoms in multilayer systems 
have been observed via mechanical means in atomic beam 
transmission and quantum reflection experiments 0, 
and via spectroscopic means y( , inter alia frequency mod- 
ulated selective reflection spectroscopy £|. 

As long as the atom-body separation is sufficiently 
large compared with the atomic radius on the one hand, 
and the typical distance between the atomic constituents 
of the body on the other hand, the vdW force can be cal- 
culated within the frame of macroscopic electrodynamics. 
A unified theory that covers both the nonretarded dis- 
tance regime — already studied by Lennard-Jones in 1932 
[f| — and the retarded distance regime was first given in 
1948 by Casimir and Polder [6j. With this in mind, the 
force has also been called Casimir-Polder force. Casimir's 
and Polder's theory is based on exact quantum electro- 
dynamics (QED), the electromagnetic field being quan- 
tized in terms of normal modes. The coupling energy of 
a ground-state atom with the body-assisted electromag- 
netic vacuum field is calculated in lowest order of pertur- 
bation theory, and the vdW force emerges as the gradient 
of this coupling energy — the vdW potential. This formal- 
ism first applied to the case of an atom placed in front 
of a perfectly conducting plate was later extended to ex- 
cited atomic energy eigenstates Q as well as to an atom 
between two perfectly conducting plates @. Moreover, 
the concept has been used to calculate the vdW force 
acting on an atom placed in front of a semi-infinite di- 
electric half spa ce M or between two dielectric plates of 
finite thickness 10]. Recently, the ideas of Casimir and 
Polder have been gen eralized to allow for dispersing and 
absorbing bodies [ill Ha ] , which inhibit electromagnetic- 
field quantization via a standard normal-mode expansion 



in general. 

In parallel with the sophistication of Casimir's and 
Polder's concept based on exact quantum electrodynam- 
ics, a semiphenomenological approach to the problem of 
the vdW force has been established and widely used. 
According to this approach, the coupling energy is ex- 
pressed in terms of correlation functions for the atom 
and/or the electromagnetic field which in turn are re- 
lated to susceptibilities via the dissipation-fluctuation 
theorem. The result — which in principle applies to ar- 
bitrary geometries — was first applied to a ground-state 
atom placed in front of a perfectly conducting half space 
[l3j . a dielectric half space (li|. and a dielectric two- 
layer system Later, atoms in excited energy eigen- 
states were included in the concept m| . Effects of surface 
roughness [17|. finite temperature 01 an d — m the case 
of the semi-infinite half space — different materials such 
as biref ring ent dielectric [19j and even magnetodielectric 
matter [2(j have been considered. 

Apart from the two main routes outlined above, a num- 
ber of other methods have been suggested and applied to 
various systems. The vdW potential of a two-level atom 
in front of a perfectly conducting half space has been de- 
rived upon using nonperturbative spectrum-summation 
techniques |2lT ] and classical electrodynamics with ran- 
dom fluctuations The problem of the vdW force 
acting on a ground-state atom in front of a dielectric 
half space has been treated via microscopic models p3j. 
S'-matrix formalism [24j, and source theory [2^. Elec- 
trostatic methods applicable to the nonretarded distance 
regime have been used to determine the vdW force act- 
ing on an excited-state atom in front of a semi-infinite 
half space filled by a birefringent dielectric [2(j, and the 
problem of the vdW force acting on an atom in front of a 
nondispcrsivc dielectric three- layer system has been stud- 
ied [23 • Within the frame of macroscopic quantum elec- 
trodynamics, a dynamical approach to the vdW force has 
recently been developed in order to study time-dependent 
forces in the case of atoms initially prepared in an arbi- 
trary excited quantum state |12| . 
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In the large body of work on vdW forces and related 
electromagnetic forces (such as the vdW force between 
two atoms or the Casimir force between two macroscopic 
bodies) the electric properties of the involved material 
objects have typically been the focus of interest. Never- 
theless, the interaction of objects possessing also notice- 
able magnetic properties — a problem which has regained 
topicality due to the recent fabrication of metamaterials 
with controllable electromagnetic properties in the mi- 
crowave regime [2^, — has been of interest. The fact 
that Maxwell's equations in the absence of (free) charges 
and currents are invariant under a duality transformation 
between electric and magnetic fields can be exploited to 
extend the notion of forces acting on electrically polariz- 
able objects to objects with magnetic properties. Thus, 
knowing the attractive vdW force between two electri- 
cally polarizable particles (e.g., atoms), one can infer 
the existence of an analogeous attractive force between 
two magnetically polarizable particles, which may be ob- 
tained from the former by replacing the electric polariz- 
abilities by the corresponding magnetic ones. In contrast, 
the force between two polarizable particles of opposed 
type is repulsive [3(j ■ While the repulsive vdW potential 
in the retarded limit obeys the same 1/r 7 power law (r, 
distance between the particles) as the attractive vdW po- 
tential (in the case of two particles of the same type), but 
is smaller in magnitude than the latter by a factor of 7/23 
[3TI l32| , the leading contribution to the repulsive vdW 
potential in the nonretarded limit is proportional to 1/r 4 , 
which contrasts with the 1 /r 6 -dependence of the attrac- 
tive vdW potential. This difference can be understood 
by regarding the first particle as an oscillating electric 
dipole creating an electromagnetic field which acts on the 
second, electrically or magnetically polarizable particle. 
Due to the fact that in the nonretarded limit the elec- 
tromagnetic field is dominated by the electrostatic field, 
the force on a electrically polarizable particle is stronger 
than the force on a magnetically polarizable one j^]. 

Similar considerations can also be made for other sys- 
tems. So, the attractive Casimir force between two in- 
finitely permeable plates corresponds to the force be- 
tween two perfectly conducting plates by virtue of du- 
ality, whereas the force between two plates of different 
type is repulsive and smaller than the equal-type force 
by a factor of 7/8 |34j[. For realistic plates the situa- 
tion becomes more involved. In particular, the repulsive 
Casimir force between a purely dielectric and a purely 
magnetic plate observed in the retarded limit shows the 
same 1/z 4 power law (z, separation of the plates) as the 
attractive force between two dielectric plates, whereas 
in the nonretarded limit the repulsive force behaves like 
1/z and the attractive force like 1/z 3 [3f|. If one of the 
plates is dielectric while the other one is magnetodielec- 
tric, then the electric and magnetic properties of the sec- 
ond plate compete in determining the sign of the Casimir 
force miS- 

It is known that the force acting on a magnetically po- 
larizable particle in front of a perfectly conducting plate 



is repulsive in the retarded limit |32j. By virtue of dual- 
ity a corresponding repulsive force is expected to act on 
an electrically polarizable particle such as a ground-state 
atom which is located in front of an infinitely permeable 
plate. Thus the question arises what kind of force could 
be observed in the case of a genuinely magnetodielec- 
tric plate. Maybe the effect of a repulsive force compo- 
nent in such a system is easier accessible to experimental 
verification than that of a repulsive component of the 
Casimir force between two macroscopic bodies, where 
force measurements are currently restricted to distance 
regimes of purely attractive forces j^. Moreover, the 
recently reported production of metamaterials with con- 
trollable magnctodielcctric properties in the microwave 
regime j^, Ei| opens the perspective of engineering vdW 
potentials with desired properties. 

In this paper we consider the vdW interaction of a 
ground-state atom with planar, dispersing, and absorb- 
ing magnetodielectric bodies. Starting from the general 
expression for the vdW potential in case of an arbitrary 
planar multilayer system, as can be derived in lowest or- 
der of perturbation theory within the frame of QED in 
linear, causal media, we give a detailed analysis of the 
vdW potential of the atom being located (i) in front of a 
magnetodielectric plate and (ii) beween two magnetodi- 
electric plates. In particular, we address the question if 
and how the competition of electric and magnetic prop- 
erties of the material can give rise to a repulsive force. 
In this context we study the influence of effects such as 
material absorption, finite layer thickness, and multiple 
reflections. 

The paper is organized as follows. In Sec. [H] the vdW 
potential of a ground-state atom in an arbitrary planar 
magnetodielectric multilayer system is given. A detailed 
analysis of typical examples is given in Sec. IIIII followed 
by a summary and concluding remarks in Sec. IIVI 

II. BASIC EQUATIONS 

Consider a neutral, nonpolar, ground-state atomic sys- 
tem such as an atom or a molecule (briefly referred to as 
atom in the following) at position ta within an arbitrary 
arrangement of linear magnetodielectric bodies, which is 
characterized by a permittivity e(r,uj) and a permeabil- 
ity /Lt(r,(x>), which are spatially varying, complex- valued 
functions of frequency, with the corresponding Kramers- 
Kronig relations being satisfied. The position-dependent 
fluctuations of the body-assisted electromagnetic field 
give rise to a force on the atom which, within leading- 
order perturbation theory, can be derived from the vdW 
potential [l2| 



U(r A ) = 




duuV 0) (m)TrG (1) (r A ,r A ,iu) (1) 



according to 

F(r A ) = -V A t/(r A ) (2) 



3 



I = 


I = 1 








/ = n 


e (w) 


£l{uj) 
























di 




dj 






z - 


= 


- * 

z 



FIG. 1: Sketch of the planar multilayer system. 
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a (0) (w) = lim 



\ " 



o 3h ^—^ w| — u> 2 — iujt 



|d 0fc | 2 



(3) 



is the ground-state polarizability of the atom in lowest 
nonvanishing order of perturbation theory [u>ko = (Ek 
— Eo)/h, (unperturbed) atomic transition frequencies; 
d fc = (0|d|fc), atomic electric-dipole transition matrix el- 
ements], and G^'(r, r',iu) is the scattering part of the 
classical Green tensor of the electromagnetic field, 



G(r,r» = G (0) (r,r» + G (1) (r,r', 



(4) 



[G'°)(r, r',u/), bulk part], which is the solution to the 
equation 



G(r,r» = 5(r-r') 



(5) 

[k(y, uf) = \i 1 (r,w)] together with the boundary condi- 
tion 



G(r, r', ijj) — > for |r — r' 



CO. 



(6) 



In what follows we assume that the bodies surrounding 
the atom form a planar multilayer system, i.e., a stack 
of n + 1 layers labelled by I (I — 0, . . . , n) of thicknesses 
di with planar parallel boundary surfaces, where e(r, u) 
— £i(u>) and /j,(t,ui) — (J-i(oj) in layer I. The coordinate 
system is chosen such that the layers are perpendicular 
to the z axis and extend from z = to z — di for I ^0,n 
and from z = to z = — oo (oo) for I = (n), cf. Fig. ^ 
The scattering part of the Green tensor at imaginary fre- 
quencies for r and r' in layer j can be given by 



(7) 



(q _L e 2 ). Here, 

/-ill) I I ■ \ l J >o \ m ) 

G w {(\,z,z ,iu) = J 



E 



i- i+ 



e+e+e- b ^ z - z,) +e-e~e b J {z - z "> 



1 



eje a r^_e 



-6j (*+«') 



-2bjdj gfcj (z+z 1 ) 



for j > 0, where 

e s x , 

(e q = qjq, q=\q\) with 



-(ige z ± 6j-eg) 



e 3 (iu)nj(iu) 



(8) 



(9) 



(10) 



are the polarization vectors for s- and p-polarized waves 
propagating in the positive/negative z-direction, rj_ and 
rj + are the generalized coefficients for reflection at the 
left /right boundary of layer j, which can be calculated 
with the aid of the recursive relations 



(H±i 

i>i±i 



! w±i 



i±l± 



I Mil 


•-£)- 






e ';±i± 


U;±i 









V fc;±i b, J ^ \ b l±1 ~T b. 



-26 !± id, ± i 7 ,P 



i±l± 



^ e-».±i*±irf ±1± 



(11) 



(12) 



(1 = 1, ... ,j for rf_, l=j, . . . ,n-l for rf + , = r£ + =0), 



£i(iu)lJLi(iu) + <j 2 



(13) 



is the imaginary part of the z-component of the wave 
vector in layer /, and finally 



DJ = l-rj_rJ + e- 



(14) 



Let the atom be situated in the otherwise empty layer 
j, i.e., £j(iu) —fij(iu) = 1 and 



(15) 



To calculate the vdW potential, we substitute Eq. Q 
together with Eq. (JSJ into Eq. QJ, thereby omitting the 
irrelevant position- independent terms. Evaluating the 
trace with the aid of the relations 

ef -ef =ef -ef = 1, (16) 
e±.e± = l, e±.e^ = -l-2(^) 2 , (17) 
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which directly follow from Eqs. ©, (|10fl arid l|13|) . we 
realize that the resulting integrand of the q-integral only 
depends on q. Thus after introducing polar coordinates 
in the g^^y-plane, we can easily perform the angular in- 
tegration, leading to 



u ^=iBl du 



(iu) 



-2bz A 



3- 



- 1 



9 9 

Ac 



-2b(dj —za) 



•j± 

D s - 



A q 
dq b 



9 9 

q c 



(18) 



Note that Eq. JSJl and thus Eq. lfTS|) also apply to the 
case j = if do is formally set equal to zero (do = 0). 

Equation ljT%|) together with Eq. © and Eqs. ifTTjl - 
(|15|l presents the vdW potential of a ground-state atom 
within a general planar magnetodielectric multilayer sys- 
tem. Note that instead of calculating the generalized re- 
flection coefficients rj ± from the permittivities and per- 
meabilities of the individual layers via Eqs. illl|) - (|13[l (as 
we shall do in this paper), it is possible to determine them 
experimentally by appropriate reflectivity measurements 
(cf., e.g., Ref. [4(|). In the case where the atom is placed 
(in free space) in front of the multilayer system (j = n) , 
Eq. (jT%|) reduces to 

U(z A ) = ^ duu 2 a^(iu) dq\e- 2b ^ (19) 
8tt 2 J a Jo b 

+ U 2 r n- 



III. SPECIFIC EXAMPLES 

Typical features of the vdW potential of an atom in the 
case of magnetodielectric multilayer systems — in partic- 
ular the competing influence of the electric and magnetic 
properties of the layers, the effect of material absorp- 
tion, the influence of finite layer thickness or multiple 
reflections — can be already illustrated by studying rela- 
tively simple systems consisting of only a few layers. 



A. Perfectly reflecting plate 

As a preliminary investigation, let us consider the ide- 
alizing case of an atom positioned in the nth (empty) 
layer in front of a perfectly reflecting (multilayer) plate, 



i.e., 



yP I - 



1. We begin with the case 



= -1, 



V 



= +1, 



(20) 



which corresponds to the limit of a perfectly conducting 
plate — ► oo, as can be seen from Eqs. (|11|) and i|12|) 



[together with Eq. I|13|) ]. Changing the integration vari- 
ables in Eq. (|19fl according to (u, q) i— ► (u,b), we obtain 
the attractive potential 



U(z A ) 



dua (0) (iu) 



dbb 2 e~ 2bZA 



4tt 2 £ Jo 

- 2 3 / dua(°\iu)e- 2 ^/ c 
l6ir 2 e z A J 



(21) 



which is exactly the result found by Casimir and Polder 
for the potential of a ground-state atom in front of a per- 
fectly conducting plate In the long-distance (i.e., re- 
tarded) limit, za> c / uj a I W A = mm ({ a -'fco|^ = lj 2, . . .})], 
the atomic polarizability (iu) may be approximately 
replaced with its static value (0) and put in front of 
the integral, leading to 



U(z A ) 



3W°>(0) 

327T 2 £ 2:| 



(22) 



In the short-distance (i.e., nonretarded) limit, z A <C c/oj a 
[uj a = max({cL>fco|& = 1,2,.. .})], we may approximately 
set e ~ 2uZA / c = 1 in Eq. (|21|l and neglect the second and 
third terms in the square brackets to recover, on recalling 
Eq. J2J), the result of Lennard- Jones p|, 



U(z A ) 



1 



4:8ire z A 



Ok 



(o|d 2 |o) 

487T£ Za 



(23) 



In contrast, if the layer facing the atom is supposed to 
be infinitely permeable, i.e., \x n -\ — > oo, Eqs. I|ll|) and 
(H2J [together with Eq. JT3J|] lead to 



<-=l, <_ = -l, (24) 
and Eq. (|19|) yields the repulsive potential 

h 



U(z A ) 



, , , , dua^(iu)e- 2uz ^ c 
lo7r £q^a Jo 



1+ 2(^) + 2(^)- 



(25) 



In particular, in the long-distance limit we have [cf. 
Eq. 



U(z A ) 



3ftco( (0 >(0) 

327T 2 £ ^A 



(26) 



which by means of a duality transformation [a^ (0) = 
ai QS> (0) t—y aif (0]j can be transformed to the result ob- 
tained in Ref. [32j for a magnetically polarizable particle 

[of polarizability aL°' (0)] in front of a perfectly conduct- 
ing plate. Application of the duality transformation to 
Eq. I|25|) generalizes the result in Ref. [32| to arbitrary 
distances. 
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B. Infinitely thick plate 



and 



To be more realistic, let us first consider an atom in 
front of a sufficiently thick magnetodielectric plate which 
may be modelled by a semi- infinite half space [n = j = 1, 
£i(u;)=/ii(a;) = l, So(w) = e(u>), ^o(co) = fj,(u>)]. Substitut- 
ing the reflection coefficients as follow from Eqs. and 
into Eq. (|T§|) . we find, on recalling Eq. |JT3J|, that (b 
= b M ) 



u(za)=! bl du 

[i(iu)b — 6m 



[j,(iu)b + b 



uV°> 



1 



M 



(iu) 



9 1 

q c 



b 



e(iu)b — b 



M 



u 2 ) e(iu)b + 6m 



(27) 



Equation l|27|l is equivalent to the result derived in 
Ref. |2(il | semiphenomenologically within the frame of lin- 
ear response theory. Note that the concept of linear re- 
sponse theory may render erroneous results when try- 
ing t o go beyond perturbation theory (cf. the remark in 
Ref. E3). 

To further evaluate Eq. (|27|l . let us model the permit- 
tivity and (paramagnetric) permeability, respectively, by 



and 



e(u) = 1 + 



iu;7 



Lu\ m - LU 2 - ilVJ n 



(28) 



(29) 



It can then be shown that in the long-distance limit, i.e., 



z A > c/lu a , z A > c/uj m [lu m -- 
reduces to (see Appendix EJ 

U(z A ) = 

where 

3feq(°)(0) 



min(o; Tc , w Tm )], Eq. 

(30) 



9± 

,4 ' 



f 



e(0)v - yje(0)n(0) - 1 + v 2 



s(0)v + i/ £ (0)ai(0) - 1 + v 2 
1 fi{0)v - y/e(0)fi(0) - 1 + v 2 



« 4 fi(0)v + y/e(0)fi(0) - 1 



(31) 



while in the short-distance limit, i.e., z A c/[w^n(0)] 
a nd/or za < c/[o; M n(0)] [w M = max(w Tc , WTm), n(0) = 
x/e(0)/i(0)], Eq. J23 leads to (see Appendix EJ 



^(za) = 3- + — , 



ZA 



where 
C 3 



167r 2 e 



du (iu 



e(iu) — 1 
e(iu) + 1 



> 



(32) 



(33) 



d = 



16tt 2 



du u 2 a^(iu) 



e(iu) — 1 



/x(iu) — 1 
fi(iu) + 1 



e(iu) + 1 
2£(iu)[e(iu)^(m) — 1] 
[e(iu) + l} 2 



>0. (34) 



We have numerically checked the asymptotic behaviour 
given by Eqs. I|30fl - H34() for the case of a two-level atom. 
From the derivation it is clear that Eqs. (|30|l and 1M2I) 
also remain valid when — in generalization of Eqs. I|28|) 
and Ij29(l . respectively — more than one matter resonance 
is taken into account. Needless to say that the minimum 
w M and the maximum w M are then defined with respect 
to all matter resonances. 

Inspection of Eq. (|31|l reveals that the coefficient C4 
in Eq. (|3"0|) for the long-distance behaviour of the vdW 
potential is negative (positive) for a purely electric (mag- 
netic) plate, corresponding to an attractive (repulsive) 
force. For a genuinely magnetodielectric plate the sit- 
uation is more involved. As the coefficient C4 mono- 
toncously decreases with increasing e(0) and monoto- 
neously increases with increasing /i(0), 



de(0) 



<0, 



>o, 



(35) 



the border between the attractive and repulsive poten- 
tial, i.e., C4 = 0, can be marked by a unique curve in 
the e(0)^(0)-plane [curves (a) in Fig. 0. In particu- 
lar, in the limits of weak and strong magnetodielectric 
properties the integral in Eq. I|31[) can be evaluated an- 
alytically. For weak magnetodielectric properties, i.e., 
Xc(0) = e(0) -1«1 and x m (0) = /z(0) -1«1, the linear 
expansions 



e(0)v - ^e(0)/i(0) - 1 + v 2 
e(0)v- 



'1 1 

2 ~ iv 2 



Xe(0) - ^Xm(O) 



(36) 



and 



fi(0)v - y/e(0)u(0) -1 + v 2 



H(0)v + y/e(0)fi(0) -l + v 2 



1 

4^2 



Xm(0) (37) 



lead to 



C 4 = 



hca^(0) 
640tt 2 £o 



[23 X c(0)-7 Xm (0)]. (38) 



For strong magnetodielectric properties, i.e., e(0) 3> 1 and 
/i(0) 3> 1, we may approximately set, on noting that large 
values of v are effectively suppressed in the integral in 
Eq. (5TJ|, 



Ve(O)MO) - l + v 2 ~ v/e(0V(0), 



(39) 
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pulsive vdW potential proportional to 1/za is predicted, 
in practice the attractive l/z^ term will always domi- 
nate for sufficiently small values of z&, because of the 
always existing electric properties of the plate. Hence 
when in the long-distance limit the potential becomes 
repulsive due to sufficiently strong magnetic properties, 
then the formation of a potential wall at intermediate 
distances becomes possible. It is evident that with de- 
creasing strength of the electric properties the maximum 
of the wall is shifted to smaller distances while increasing 
in height. 

In the limiting case of weak electric properties, i.e., 
cjpc/tJTc^l and wp c /wp m <l [recall Eqs. l(2"5)l and lt29l) ] 
one can thus expect that the wall is situated within the 
short-distance range, so that Eqs. 1(32(1 - 1(34(1 can be used 
to determine both its position and height. From Eq. i(32|) 
we find that the wall maximum is at 



max 
Z A 



3C 3 



(41) 



and has a height of 



FIG. 2: Border between attractive and repulsive long- 
distance vdW potentials of an atom in front of (a) a thick 
and (b) a thin magnetodielectric plate according to Eqs. 13 H 
(C4 = 0) and 16UH (D 5 — 0). The broken curves show the 
asymptotic behaviour as given by Eqs. £g) (inset) and JIOJl 
in case (a) and by Eqs. 1611 (inset) and 16211 in case (b). 



thus 

C 4 = 



64tt 2 £ 
1 4 

Z ~ 3 ~ 



-z)- 



Z + 2Z 1 - 2Z i \n \ 1 



2 4 

Z 2 " + ~Z 
1 

z 



ln(l + Z) 



(40) 



where Z = y/ /i(0)/e(0) is the static impedance of the 
material. Setting 64 = in Eqs. 1(38(1 and 1(4011 . we obtain 
the asymptotic behaviour of the border curve in the two 
limiting cases. In particular, from Eq. I|40|) it follows that 
Z = 2.26. In conclusion one can say that in the long- 
distance limit a repulsive vdW potential can be realized 
if the static magnetic properties are stronger than the 
static electric properties, Xm(0)/xe(0) > 23/7 = 3.29 for 
weak magnetodielectric properties, and /i(0)/e(0) > 5.11 
for strong magnetodielectric properties. 

Apart from the different distance laws, the short- 
distance vdW potential, Eq. 1)32(1 . differs from the long- 
distance potential, Eq. (|30|l . in two respects. First, the 
relevant coefficients C3 and C\ are not only determined 
by the static values of the permittivity and the perme- 
ability, as is seen from Eqs. I|33|l and (|34|l . and second, 
Eqs. J22J-(|23J reveal that electric and magnetic proper- 
ties give rise to potentials with different distance laws 
and signs [C3 > dominant (and C\ > 0) if e ^ 1 and \x = 
1, while C3 = and C\ > if e = 1 and fi ^ 1]. However, 
although for the case of a purely magnetic plate a re- 




(42) 



In order to estimate the integrals in Eqs. 1(33(1 and 1341) for 

the coefficients C3 and C\ , respectively, let us restrict our 
attention to the case of a two-level atom and disregard 
absorption (7 C ~ 0, 7 m ~ 0). Straightforward calculation 
then yields (wp c /wt c -C 1, wp c /wp m <C 1) 
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967T£ Wf c Wio + ^Tc 
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and 



16tt 2 



duu 2 a^ (iu) 



n(iu) — l /j,(iu) — l 
[i(iu) + l 2 



^ol d Ol| 2 ^Pm ^10( 2 ^10 +^Sm +^Tm) 
967T (wio + WSm) (Wio + WTm) 



(44) 
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(|44|l into Eqs. I(41|) and 142(1. respectively, eventually leads 



1, ,2 



1/21 



Substitution of Eqs. I|43() and 



to 



WPm WTc y Wlo(Wl0 + U)Te) 
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Y (2cJio +WSm +W Tm ) 



(45) 
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FIG. 3: The vdW potential of a ground-state two-level atom 
situated in front of an infinitely thick magnetodielectric plate 
is shown as a function of the distance between the atom and 
the plate for different values of n(0) (ujp c /ujio — 0.75, u>re/wio 
= 1.03, WTm/wio = 1, Jc/loio =7m/^10 =0.001). 



FIG. 4: The vdW potential of a ground-state two-level atom 
situated in front of an infinitely thick magnetodielectric plate 
is shown as a function of the distance between the atom and 
the plate for different values of wm [m(0) = 5, u}p c /uiio = 0.75, 
WTe/wio = 1.03, 7c/^io =7m/wio =0.001]. 



Note that consistency with the assumption of the wall 
being observed at short distances requires that z™ ax -C 
c/u/jyj — a condition which can be easily fulfilled for suffi- 
ciently small values of wp /wp m - Inspection of Eq. H46|) 
shows that the height of the wall increases with increas- 
ing wp m , decreasing WT m , and decreasing wp G /o;Tc = 
\/e(0) — 1. Since the dependence of C/(z™ ax ) on ujp m 
is much stronger than its dependendence on WTm, the 
wall height increases with <x>T m for given wp m /wTm = 

The distance dependence of the vdW potential, as cal- 
culated from Eq. (|27() for a two-level atom in front of 
a thick magnetodielectric plate whose permittivity and 
permeability are modelled by Eqs. (|2"S)l and respec- 
tively, is illustrated in Figs. |3] and 0] The figures re- 
veal that the results derived above for the case where 
the potential wall is observed in the short-distance range 
remain qualitatively valid also for larger distances. So, 
from Fig. |3 it is seen that, for chosen values of WTm and 
7 m , the potential wall begins to form and grows in height 
as /z(0) increases, while Fig. 0] confirms that, for chosen 
values of /i(0) and j m , the height of the wall increases 
with WTm- In conclusion one can say that the formation 
of a noticeable potential wall requires materials whose 
static permeability substantially exceeds the static per- 
mittivity, thereby featuring magnetic resonance frequen- 
cies as high as possible. 

To study the dependence of the vdW potential on ma- 
terial absorption as characterized by the parameters 7 C 
and 7 m in Eqs. 128|) and (|29|l . we first consider the lim- 
iting behaviour of the potential for long and short dis- 
tances. As the potential in the long-distance limit can be 
given in terms of the static permittivity and permeabil- 



ity, which do not depend on the absorption parameters, 
material absorption has no influence on the vdW force 
for asymptotically large distances. In contrast, absorp- 
tion can affect the potential in the short-distance limit. 
From Eqs. l(2"8|) and l(2"9")l the inequalities 



de(iu) 



<0, 



dfi(iu) 

37m 



< 



(47) 



are seen to be valid. Combining them with Eqs. i.VMi and 
(|34|l reveals that 



dC 3 


<o, 


dC 3 




97m 


dCi 


<o, 


dd 


97c 





= o, 

< 0. 



(48) 
(49) 



Provided that the magnetic properties of the medium 
are sufficiently strong to support the formation of a po- 
tential wall, these inequalities imply [cf. Eq. (|32|l ] that 
increasing j c (7 m ) leads to a shift of the wall towards 
smaller (larger) distances, while increasing (decreasing) 
its height. Thus an increase of j e yields a stronger repul- 
sive potential, whereas a simultaneous increase of both 
absorption parameters is expected to lead to a reduction 
of the wall height in general. This behaviour is confirmed 
by the examples shown in Fig. [31 where the vdW poten- 
tial of a two-level atom as given by Eq. (|27|l is displayed 
as a function of the distance between the atom and the 
plate for different values of the two absorption parame- 
ters. Note the reduced influence of absorption at large 
distances — in agreement with the arguments given above. 



In view of left-handed materials (see, e.g., Refs. [28, 29, 
l4l|). which simultaneously exhibit negative real parts of 
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[^ 2 (m)6 2 — 6^] tanh(6M^) 
2[i(iu)bb M + [fJ. 2 (iu)b 2 + b^] tanh(6 M d) 



(50) 



It is obvious that the integration in Eq. (|50|) is effectively 
limited by the exponential factor e~ 2fcZA to a circular re- 
gion where 6<1/(2za). In particular, in the limit of a 
sufficiently thick plate, d^$> z A , the estimate 



FIG. 5: The vdW potential of a ground-state two-level atom 
situated in front of a infinitely thick magnetodielectric plate 
is shown as a function of the distance between the atom and 
the plate for different values of the absorption parameters: (a) 
7eMo = 0.001, j m /u)io = 0.001, (b) 7c /o; 1 o = 0.001, 7m /aji 
= 0.05, (c) 7e /wio=0.05, 7 m /wio = 0.001, (d) 7c /oj 10 = 0.05, 
7m/i^io=0.05 (u>p c /uj 10 =0.75, wtc/^io = 1-03, wp m /wio=2, 
WTm/wio = 1). 



e(w) and within some (real) frequency interval such 
that the real part of the refractive index becomes nega- 
tive therein, the question may arise whether these ma- 
terials would have an exceptional effect on the ground- 
state vdW force. The answer is obviously no, because the 
ground-state vdW potential as given by Eq. (|27|l is ex- 
pressed in terms of the always positive values of the per- 
mittivity and the permeability at imaginary frequencies. 
Clearly, the situation may change for an atom prepared 
in an excited state. In such a case, the vdW potential is 
essentially determined by the real part of the Green ten- 
sor taken at frequencies close to the transition frequencies 
of the atom [lj| ■ When there are transition frequencies 
that lie in frequency intervals where the material behaves 
left-handed, then particularities may occur. 



by^d > bd 



Plate of finite thickness 



Let us now consider an atom in front of a magnetodi- 
electric plate of finite thickness d\ = d [n = j = 2, Si(uj) 
= s(u>), Mi(w) = £ (oj) = e 2 (w) = 1, Mo(w) = M2M 
= 1]. Substituting the reflection coefficients calculated 
from Eqs. (fTT|) and (|T^|l into Eq. ljl"§|) . we derive (fei = 



2z A 



> 1 



(51) 



[recall Eqs. I|13l) and (|15fl ] is approximately valid within 
(the major part of) the effective region of integration, 
and one may hence make the approximation tanh(6M^) 
~ 1 in Eq. (|5"U|) . which obviously leads back to Eq. lj2*7|) 
valid for an infinitely thick plate. On the contrary, in the 
limit of an asymptotically thin plate, n(0)d<^z A , we find 
that the inequalities 



by[d < \/e(iu)^i(iu) bd < \f e(0)^(0) bd 
n(0)d 



< 



2z A 



< 1 



(52) 



hold in the effective region of integration, and one may 
hence perform a linear expansion of the integrand in 
Eq. I|50(l in terms of ^m^, resulting in 



U(z A ) 



hfigd 



8tt 2 
- (1 + 2 



du u 2 a^ (iu) 
Jo 



dq 1 e- 2bZA 
b 



q 2 c 2 \ s 2 (iu)b 2 — 6^ ^ 2 {iu)b 2 — b 2 



M 



2e{iu)b 



2fi(iu)b 



(53) 



Provided that the magnetic properties are sufficiently 
strong, the formation of a repulsive potential wall can be 
also observed in the case of a genuinely magnetodielectric 
plate of finite thickness. Typical examples of the vdW 
potential obtained by numerical evaluation of Eq. I|50l) 
for a two-level atom are shown in Fig. In the figure, 
the medium parameters correspond to those which have 
been found in Sec. IIII Bl to support the formation of a 
potential wall in the case of an infinitely thick plate. We 
see that the qualitative behaviour of the vdW potential 
is independent of the plate thickness. In particular, all 
curves in Fig. [B] feature a repulsive long-range potential 
that leads to a potential wall of finite height, the poten- 
tial becoming attractive at very short distances. How- 
ever, the position and height of the wall are seen to vary 
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FIG. 6: The vdW energy of a ground-state two- level atom 
situated in front of a magnetodielectric plate is shown as a 
function of the distance between the body and the interface 
for different values of the plate thickness d {up c /u\o =0.75, 

WTe/wio = 1.03, Wp m /Wl0 = 2, LJTm/^10 = 1, Je/UlO = Jm/UxO 
= 0.001). 



with the thickness of the plate. While the position of the 
wall shifts only slightly as the plate thickness is changed 
from very small to very large values, the height of the 
wall reacts very sensitively as the plate thickness is var- 
ied. For small values of the thickness the potential height 
is very small, it increases towards a maximum, and then 
decreases asymptotically towards the value found for the 
infinitely thick plate as the thickness is increased further 
towards very large values. It is worth noting that there 
is an optimal plate thickness for creating a maximum 
potential wall. In this case the magnitude of the plate 
thickness is comparable to the position of the potential 
maximum — a case which is realized between the two ex- 
tremes of infinitely thick and infinitely thin plates. 

In order to gain further insight into the competing elec- 
tric and magnetic effects on the formation of an poten- 
tial wall, let us study the case of an asymptotically thin 
plate as described by Eq. I|53|) in more detail and com- 
pare it with the case of an infinitely thick plate studied in 



Sec. IIII Bl In the long-distance limit, za^> c/uj a , c/uu 
Eq. iftnfy reduces to (see Appendix El 



Z A 



where 
Da = 



hca^(0)d 
160tt 2 £ 



14e 2 (0) - 9 6^ 2 (0) - 1 



e(0) 



M0) 



(54) 



(55) 



while in the short-distance limit, za <Cc/[o^n(0)] and/or 
Za ^ c /[ w M n (0)L -^q. can be approximated by (see 



Appendix [XJ 
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64tt 2 £ Jo 



and 



64tt 2 



du u 2 a.(°\iu)- 



e(iu) 



e 2 {iu) - 1 
e(iu) 



(56) 



(57) 



[i 2 (iu) — 1 2[e(iu)fj,(iu) — 1] 



e(iu) 



>0. (58) 



In the case of an asymptotically thin plate the border 
between attractive and repulsive potentials is determined 
by the equation D§ = 0, because Eq. (|55|l reveals that 



0D 5 
de(0) 



<0, 



dD 5 
Bn{0) 



> 



(59) 



[cf. Eq. I|35|) valid for an infinitely thick plate]. Since 
for an asymptotically thin plate — in contrast to the in- 
finitely thick plate — the influence of electric and mag- 
netic properties can be completely separated into a sum 
of two terms, the equation D§ — can be solved analyti- 
cally, leading to 



M(0) 



14e 2 (0) - 9 + \/l96e 4 (0) - 228e 2 (0) + 81 
12e(0) 



(60) 



[curves (b) in Fig. For sufficiently weak magne- 

todielectric properties, i.e., Xe(0) =e(0) — 1 <C 1, x m (0) = 
/i(0) — 1 <C 1, a linear expansion of the right-hand side of 
Eq. (|60|l reveals that a repulsive vdW potential can be re- 
alized if the static magnetic properties are stronger than 
the static electric properties by a factor Xm(0)/Xe(0) > 
23/7 = 3.29, corresponding to 



D 5 



hca^(0)d 
160tt 2 £ 



[23 Xe (0)-7xm(0)], (61) 



as can be seen by linearly expanding the right-hand side 
of Eq. (|55|l . By comparing Eqs. (|3*8*|l and lfHT)l we real- 
ize that in the limit of weak magnetodielectric properties 
the border between attractive and the repulsive vdW po- 
tentials is the same for the infinitely thick plate and the 
asymptotically thin plate (cf. the inset in Fig. [21 . This 
result is an immediate consequence of the fact that in 
this case the thick-plate potential is a linear superposi- 
tion of thin-plate potentials [see Sec. IIII Dl Eq. (|72|l ]. For 
strong magnetodielectric properties, e(0) ^> 1, fi(0) ^> 
1, the asymptotic behaviour of the right-hand side of 
Eq. (|60|l shows that the vdW potential becomes repul- 
sive if /z(0)/e(0) > 7/3 = 2.33, corresponding to 



W°)(0)d r , s , o 
A = -^^[7 £ (0)-3 M (0)], 



(62) 
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which follows from the corresponding asymptotic expan- 
sion of the right-hand side of Eq. i|55|). Hence the re- 
gion in the e(0)^(0)-plane that corresponds to a repul- 
sive vdW force is slightly increased in comparison to the 
infinitely thick plate. 

As in the case of an infinitely thick plate, the elec- 
tric and magnetic properties of the medium give rise 
to competing attractive and repulsive potential compo- 
nents, where again the attractive potential component re- 
sulting from the electric properties dominates in the limit 
za [see Eqs. (|56|) - (I58|) ]. This implies that for weak 
electric properties (wp e /wTe *C 1, wp c /wp m «1) a po- 
tential wall is formed in the short-distance range. From 
Eq. ffify it then follows that the wall is situated at 



max 

Z A — 



and has a height of 



u(zT 



2D 4 

~d7 



Dl 
4L> 4 



(63) 



(64) 



For a two-level atom interacting with an almost nonab- 
sorbing (7 e ~0, 7 m ~0) single-resonance medium exhibit- 
ing weak electric properties (wpe/wTe^l, wp c /wp m <l), 
the coefficients D4, Eq. I|57|l . and D2, Eq. I|58|) . can be 
evaluated according to 



Comparing Eqs. JH^I and (JHSJ) with Eqs. (03 and (gBJ) 

valid for an infinitely thick plate, we find that the depen- 
dence of both the position and the height of the potential 
wall on the material parameters is very similar, so that 
the criteria for having a noticeable potential wall given 
below Eq. Ij46|l also apply to the case of an asymptot- 
ically thin plate. From the result that the position of 
the wall is almost the same in both cases it may be ex- 
pected that the wall position slowly varies with the plate 
thickness, which is in full agreement with the exact re- 
sults in Fig. [fj] Further, the height of the wall is — in 
agreement with Fig. [5| — considerably smaller for the as- 
myptotically thin plate. This can be seen by applying 
d/z% ayL < ^/e(0H0)d/z£ lax < 1 together with Eq. in 
Eq. leading to 



u(zT 



101 



UJp e 



7687T£ C 3 

wio(4wi 



^Tc 
3wLm + 



3(^10 + W Lm )(wi +WTm) 



(69) 



the right-hand side of which is comparable to the right- 
hand side of Eq. I|46l) . Recall that the wall height does not 
monotonously increase with the plate thickness in general 
[see Fig.|S], as could be expected from comparing the two 
limiting cases. 
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96tt 
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(66) 



(uj Lm = sj uo\ m + u>p m ) ■ Substitution of Eqs. JHSJ) and JHEJ 
into Eqs. (J^J and leads to 



max 
Z A 




(67) 



(with the consistency requirement z™ ax <C c/^m being 
fulfilled for sufficiently small values of cjp G /wp m ) and 



u( Z r x ) 



dldoi 



J Pm W Te W 10 + w Te 



11527T£ C 4 
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-1 2 
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(68) 



D. Power laws and medium-assisted correlations 

Comparing the asymptotic power laws l|30|l and 1321) 
obtained for an infinitely thick plate with those obtained 
for an asymptotically thin plate, Eqs. lfS^|) and we 
see that in the latter case the powers of 1/za are univer- 
sally increased by one. In both cases the long-distance 
vdW potential follows a power law that is independent of 
the material properties of the plate, the sign being deter- 
mined by the relative strengths of magnetic and electric 
properties (a purely electric plate creates an attractive 
vdW potential, while a purely magnetic plate gives rise 
to a repulsive one). Further, the short-distance results 
for plates of different material properties differ in both 
sign and leading power law (the repulsive potential cre- 
ated by a purely magnetic plate being weaker than the 
attractive potential created by a purely electric plate by 
two powers in the atom-plate separation). It is worth 
noting that a similar behaviour, i.e., the same hierarchy 
of power laws and the same signs have been found for 
the vdW force between two atoms (^J 132I l33l| and for 
the Casimir force between two semi-infinite half-spaces 
[35j]. This is illustrated in Tab. HJ where the asymptotic 
power laws found for an atom interacting with an in- 
finitely thick plate, Eqs. 130fl and H32|) . and an asymp- 
totically thin plate, Eqs. and (JSIjJ, are summarized 
and compared to those obtainable for the interactions 
between two atoms or two half-spaces, respectively. 

For weak magnetodielectric properties, i.e., Xc(* u ) = 
e(iu) — 1 1 and Xm(iu) = fi(iu) — 1 -C 1, the similarity 
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number of thin- plate vdW potentials, 



TABLE I: Signs and asymptotic power laws of the forces 
between various polarizable objects. In the table heading, e 
stands for a purely electric object, m denotes a purely mag- 
netic one. The signs + and — denote repulsive and attractive 
forces, respectively. 
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A 1 U(z A ) = / — A 1 U d (z). (72) 
Jz A d 



For media with stronger magnetodielectric properties, 
medium- assisted correlations prevent vdW-type forces 
from being additive. This can be demonstrated by ex- 
panding the vdW potentials given by Eqs. 127(1 and 1531) 
to second order in % e (m) and Xm(iu), resulting in 
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of the results shown in Tab. [I] can be regarded as being a 
consequence of the additivity of vdW-type interactions. 
In fact, in this case (which for a gaseous medium of given 
atomic species corresponds to a sufficiently dilute gas) all 
results of the table can be derived from the vdW inter- 
action of two single atoms via pairwise summation. The 
additivity can explicitly be seen when comparing the re- 
sult found for the asymptotically thin plate with that of 
the infinitely thick plate. Expanding the vdW potential 
of an infinitely thick plate, Eq. (|27|l . in powers of \ e {iu) 
and x m (ra), we find that the leading, first-order contri- 
bution is given by 
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Xm(iu) 



(70) 



while the first-order contribution to the vdW potential 
of an asymptotically thin plate, Eq. H53[) reads 
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(71) 



Comparison of Eqs. (|70|) and (|71|l shows that to leading 
order in Xeiiu) and Xm(iu) the vdW potential of an in- 
finitely thick plate is simply the integral over an infinite 
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(74) 



respectively. The leading (second-order) correction due 
to medium-assisted correlations can be obtained from the 
vdW potential of an atom in front of two asymptotically 
thin plates. Physically, it can be ascribed to the process 
of radiation being reflected at the back (left) plate while 
aquiring finite phase shifts upon transmission trough the 
front (right) plate. The calculation yields (see Appendix 

03 
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Xc(iu)xm(iu) >. (75) 



Note that the leading correction due to multiple reflec- 
tions between the plates and fractional transparency of 
the front (right) plate are of third order in Xe(iu) and 
X m (iu), and are thus not relevant for the second-order 
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correction considered here. Comparing Eqs. JTSJi CD; 
and J7SJ, one can easily verify that 



A 2 U{z A ) 



CO J 

dz 
d 



A 2 U dd (z A ,s). (76) 



As a consequence of medium- assisted correlations the co- 
efficients of the asymptotic power laws in Tab. [I] cannot 
be related via simple additivity arguments in general. 
However, we note from Tab. [I] that the corrections only 
change the coefficients of the asmptotic power laws, not 
the power laws themselves. 



E. Atom between two infinitely thick plates 
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In Sees. IIII Bl and IIII CI we have seen that for suffi- 
ciently strong magnetic properties a single magnetodi- 
electric plate can feature a potential wall. This suggests 
that two such plates can feature a potential well, where 
the effect of multiple reflections between the plates must 
be taken into account. Let us consider the simplest case 
of an atom placed between two identical infinitely thick 
magnetodielectric plates which are separated by a dis- 
tance di = s [n = 2, j=l, £i(u>) = jUl(w) = 1, £q(u)) =£2(&>) 
= e(ui), ^o(w) =^2(w) = fi(u>)]. From Eq. (fTSf> together 
with Eqs. l(TT|) - ((T4l it then follows that (6 = 6 2 = 6m) 
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u 2 J D\ e(iu)b + 6m 



where the coefficients 



D{ = 1 



n(iu)b — 6m 



/j.(iu)b + b M 
e(iu)b — 6m 



e(iu)b + 6m 



e' 2bs < 1, 



e- 2bs < 1 



(77) 

(78) 
(79) 



describe the effect of multiple reflections of radiation be- 
tween the two plates, as can be seen from the expansion 



i+ c 



D -2bs 



^ — bsn 

e r. 



(80) 

As a consequence of multiple reflections the vdW poten- 
tial of an atom between the two plates, Eq. (|77|l . can 
be different from the sum of two single-plate potentials, 
Eq. J23). 

Examples of the vdW potential for a two-level atom 
between two identical infinitely thick magnetodieletric 



FIG. 7: The vdW potential of a ground-state two-level 
atom situated between two infinitely thick (a) magnetodi- 
electric plates (wpc/wio = 0.75, lutc/loio = 1.03, o;p m /wio = 2, 
cjTm/^io = 1, 7e/^io = 7m/wio = 0.001) (b) dielectric plates 
[u(uj) = 1, other parameters as in (a)], (c) magnetic plates 
[e(ui) = 1, other parameters as in (a)], which are separated by 
a distance s— 15c/ljio, is shown as a function of the position 
of the atom. 



plates as given by Eq. I|77|l are plotted in Figs. an d 
|H1 In the case of the parameters chosen in Fig. multiple 
reflections are negligible, so that the potentials effectively 
reduce to sums of single-plate potentials. This obviously 
results from the smallness of the relevant reflection coeffi- 
cients together with the relatively large distance between 
the plates. From Eqs. I|llf) and (|12|l one can easily verify 
that 

r{ ± (iu, q) < lim r s 1± (0, q) = , (81) 

ri ± (iu,q) < lim r? ± (0, ? ) = £ -^—\ . (82) 

Hence in the case of the parameters in Fig. Eta) we have 
r i- r i+ — 0.67, r\_r\ + < 0.044. In order to demonstrate 
the effect of multiple reflections, in Fig. |§]we have (arti- 
ficially) increased the reflection coeffiecients, so that al- 
most perfect reflection is realized (r*_r* + < 1 — 1 x 10 



10 



p p 



< 1 — 7.5 x 10~ 10 ), and reduced the plate separa- 
tion. It is seen that multiple reflections lead to a slight 
lowering of the vdW potential in the region near the mid- 
point between the two plates. 



IV. SUMMARY AND CONCLUSIONS 

We have studied the problem of the van der Waals 
force acting on a ground-state atom in the presence of 
planar, dispersing, and absorbing magnetodielectric bod- 
ies. Considering an arbitrary planar multilayer system 
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0.00075, r 




FIG. 8: The vdW energy of a ground-state two- level atom 
situated between two infinitely thick magnetodielectric plates 
(lu Pc /lo 10 = 0.75 x 10 5 , WTe/wio = 1-03, i* m /uio = 2 x 10 5 , 
WTm/wio = 1, 7c/<^io = 7m/wio = 0.001), which are separated 
by a distance s — Qc/ljio, is shown as a function of the position 
of the atom, Eq. (I77t . For comparison, the sum of two single- 
plate potentials according to Eq. 12711 is also displayed (dashed 
lines). 



and restricting our attention to the lowest (nonvanish- 
ing) order of perturbation theory, we have given a gen- 
eral expression for the vdW potential. The effect of the 
multilayer system is expressed in terms of generalized re- 
flection coefficients, which on their part are determined, 
inter alia, by the (complex) permittivities and permeabil- 
ities of the layers. Applying the formula to the cases of 
an atom being in front of a magnetodielectric plate and 
between two such plates, we have placed special emphasis 
on the competing attractive and repulsive force compo- 
nents associated with the electric and magnetic matter 
properties, respectively. Both numerical and analytical 
results are given, the latter referring to some limiting 
cases such as the asymptotic behaviour of the potential 
for thick and thin plates in the long- and short-distance 
limits. 

In contrast to the well-known attractive vdW force 
generated by a purely dielectric plate, a purely mag- 
netic plate leads to a repulsive force. In the case of gen- 
uinely magnetodielectric material, the influence of the 
magnetic properties can thus considerably reduce the 
strength of the vdW force and — for sufficiently strong 
magnetic properties — even create a repulsive potential 
wall of finite height. The numerical results show that the 
height of such a potential wall sensitively depends not 
only on the relative strengths of the electric and mag- 
netic properties, but also on the thickness of the plate. 
In particular, they suggest that the maximum height is 
realized in the case when the thickness of the plate is 
comparable to the distance of the potential maximum 



from the plate. Comparing the results obtained for an 
infinitely thick plate with those found for an asymptoti- 
cally thin plate, we have found striking similarities which 
for weakly magnetodielectric media can be explained by 
the additivity of vdW potentials. Moreover, we have ex- 
plicitly demonstrated how medium-assisted correlations 
lead to a breakdown of additivity for media with stronger 
magnetodielectric properties. For an atom being situated 
between two magnetodielectric plates each of which fea- 
tures a potential barrier, a potential well of finite depth 
can be formed. If the plates possess a sufficiently high 
reflectivity while being relatively close together multiple 
reflections can prevent the resulting potential from being 
simply the additive superposition of the two single-plate 
potentials. 

The results show that the advent of artificially made 
materials with controllable magnetodielectric properties 
will offer novel possibilities of realizing vdW potentials 
on demand. The provided analysis of typical effects rel- 
evant for controlling the vdW force in the case of one 
and two magnetodielectric plates — namely the compe- 
tition between electric and magnetic properties of the 
material in the formation of the potential, material ab- 
sorption, plate thickness, and multiple reflections — can 
of course be extended to more complex multilayer sys- 
tems by further evaluating the general formula for the 
vdW potential of a ground-state atom in planar multi- 
layer systems. 
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APPENDIX A: LONG- AND SHORT-DISTANCE 
LIMITS 

The long-distance (short-distance) limit corresponds to 
separations za between the atom and the multilayer sys- 
tem which are much greater (smaller) than the wave- 
lcnghts corresponding to typical frequencies of the atom 
and the multilayer system. To obtain approximate re- 
sults for the two limiting cases, let us analyze the u- 
integrals in Eqs. (|27|l and (|53|) in a little more detail and 
begin with the long-distance limit, i.e., 

za > — , z A > — , (Al) 
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where w A = min({ajfco|^ — 1; 2 . . .}) is the lowest atomic 
transition frequency, and tu^ = min(;jTe) ^Tm) is the low- 
est medium resonance frequency. For convenience, we 
introduce the new integration variable v = cb/u and trans- 
form the integral according to 



du I dq — e 



dv I du- e - 2zAVU / c . 



o 

where 6m has to be replaced according to 
u 



5M 



\J e{iu) /j,(iu) — 1 + v 2 



(A2) 



(A3) 



Inspection of Eqs. and J53J| together with Eq. (|X2|) 
reveals that the frequency interval giving the main con- 
tribution to the respective it-integral is determined by a 
set of effective cutoff functions, namely 



/(«) 



_ g-22 A «/c 



(A4) 



Combining Eqs. l|A2|l . l|Al)l . and {Hp with Eq. fgfi 
and Eq. I|53l) . repectively, and evaluating the remaining 
it-integrals we arrive at Eq. (|3U[) [together with Eq. i|31[l ] 
and Eq. (JSlJl [together with Eq. 

The short-distance limit, on the contrary, is defined by 



ZA < 



w+n(0) 



and/or za <C — ^ 



win(O) 



(A16) 



where o;^ = max({uJko\k — 1,2,.. .}) is the highest in- 
neratomic transition frequency, w M = max(wTc,WT m ) is 
the highest medium resonance frequency, and n(0) = 
y/e(0)n(0) is the static refractive index of the medium. 
Again, it is convenient to change the integration variables 
in Eqs. (|27[l and l|53[) . but now we transform according 
to 



du I dq — e 



du / dbe- 2bZA 

(i Ju/c 



(A17) 



9k{u) 



(A5) 



1 + (u/uj k o) 2 ' 

which enter via the atomic polarizability, cf . Eq. @ , and 

' (A6) 
(A7) 



h c {u) - 
h m (u) 



1 + (u/lotc) 2 ' 
1 

1 + (tl/wTm) 2 



which enter via e(iu) and n(iu), cf. Eqs. and 
The cutoff functions obviously give their main contribu- 
tions in regions, where 



< 



< 



c 
2^ 


for 


/(«), 


(A8) 




for 




(A9) 


^Tc 


for 


/l c (u), 


(A10) 


WTm 


for 


h m (u). 


(All) 



Combining Eq. (|A8|> with Eq. (|A1|> . we find that the func- 
tion /(it) effectively limits the it-integration to a region 
where 



u u . c 

< — < 7, « !> 

WfcO oj a 2za^a 

U it . c 

— < — < - — - « 1, 

utc uj m 2z A w M 
U ^ It c 

< — < < 1. 



M 



2^AW M 



(A12) 
(A13) 
(A14) 



Performing a leading-order expansion of the integrands 
in Eqs. 127(1 and (|53[1 in terms of the small quantities 
it/wfcQj u/u>t c , and u/uiT m , we may set 



a (0) (m) ~ a w (0), e(m) ~ e(0), ~ fi(0). (A15) 



v(0) 



where 6m has to be replaced according to 



— [e(iu)/j,(iu) - l] + 6 2 



(A18) 



Combining Eqs. (|A9 fl — (|Alip with Eq. (|A16|) reveals that 
the functions gk(u), h e (u), and /i m (it) limit the it- 
integration to a region where 



ZA"Uy/e(iu)n(iu) — 1 za^h(O) 



< 1 



and / or 

z A u^Je{iu)n(iu) - 1 z A ^ M n (°) 



< 1. 



(A19) 



(A20) 



A valid approximation to the it- integrals in Eqs. 127|) 
and l|53l) can hence be obtained by performing a 
Taylor exansion in ZAUy/e(iu)fi(iu) — 1/c. To that 
end, we apply the transformation l|A17(l to Eqs. I|27|) 
and (|53|l . respectively, retain only the leading-order 
terms in u\J ' e{iu)n(iu) — l/(c6) (which after carrying 
out the 6-integral will yield the leading-order terms in 
Zpju^f £(iu)n{iu) — 1/c) and obtain 



U(z A ) 



'8tt 2 



duuV )(iit) / d6e- 2b2A 

Ju/c 



x 2- 



c 2 6 2 e(m) — 1 fe(Mi) — 1 /x(m) — 1 



u 2 e(iu) + 1 [e(m) + l /i(iit) + 1 



+ 



2e(iit)[£(iii)/i(iit) — 1] 
[e(zit) + l] 2 



(A21) 
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and 



U{z A ) 



hfiod 



duu 2 a^ (iu) 



dbbe 



-2bz A 



)C 2 b 2 e 2 (m) - 1 
u 2 2e(iu) 



/c 

e 2 {iu) — 1 
2e(iu) 



H 2 (iu) — 1 e(iu)[i(iu) — 1 
2/i(iu) 



e(iu) 



(A22) 



After evaluating the 6-integrals and keeping only the 
leading-order terms in uza/c [note that Eqs. ljA9(l - l|All|) 
together with Eq. I|A16|I imply uza/c<^1], Eq. I|A21[I and 
Eq. (|A22I) . respectively, result in Eq. (|3*2)l [together with 
Eqs. (HH and ^] and Eq. ^ [together with Eqs. JS7JI 
and ijSHJl]. 



APPENDIX B: DERIVATION OF EQ.^U 

In order to derive Eq. (|75|l . we consider the vdW po- 
tential of two plates of thickness d\ = d, d 3 = d 1 which 



are separated by a distance e?2 = s [n 
e(u),e 3 {uj)=e'(uj), = /i(cj), fj, 3 (uj) : 

e 2 (^) = £4(w) = 1, /io(^) = A*2(^) = M4(w) = 1]. We as- 
sume both plates to be asymptotically thin, ^/e(0)/i(0)(i 
<C za, \/e'(0)n'(0)d' ^ Z A, so that the inequalities 6m d 
< 1, 6 M d' < 1 Oi = &m, ^3 = b' M ) are valid, cf. Eq. Q . 
Use of Eqs. {TTJ and {HI for / = ?i = 4 and Z = 3, followed 
by a linear expansion in terms of b' M d' , yields 



2/x'(m)6 " * " 2 



x 1- 



M ' 2 (m)6 2 + 6 M 
H'(iu)b 
^ 2 ^u)b 2 ~b' 2 ^ 2bSrS 



2fi'(iu)b 
e' 2 (iu)b 2 - b»;, 
2e'(iu)b 



(Bl) 



M d' + e- 2hs r? d' 



e' 2 (m)6 2 + 6 M 
e'(iu)b 
s' 2 ^u)b 2 -b' 2 ^ 2bs ^ p 



2e'(iu)b 



-e rS 



(B2) 



while use of the same equations for I = 2 and I — 1 together 
~q_ = leads to, upon linearly expanding in 



with r. 



o- 



terms of 6m d, 



fi?(iu)b 2 - 6 M d 

2fj,(iu)b 
e 2 (iu)b 2 -b 2 M d 

2e(iu)b 



(B3) 
(B4) 



Substituting Eqs. JB3) and JB4|) into Eqs. CED and JB2), 

respectively, and neglecting terms which are quadratic in 
&M<i, we arrive at 



V ,2 (iu)b 2 



2fi'{iu)b 
x 1 - 



y' 2 (iu)b 2 + b' 2 
fi'(iu)b 



2/i(m)6 



(B5) 



£ ' 2 M6 2 -6 M d , £ 2 M6 2 - 6 2 1 
2e'(iu)b 2e(iu)b 
x e' 2 (iu)b 2 + b$ d ; 
e'(iu)b 



- 2bs d 



(B6) 



The leading correction due to medium correlations can 
be extracted from Eqs. 1|B5(1 and (|B6|) by retaining only 
the two-plate contribution, i.e., the term which is lin- 
ear in both bud and b' M d' . We expand the result up to 
linear order in Xe(iu), Xm(iu), x' c {^ u ) =s'{iu) — 1, and 
x' m (iu) = fj,'(iu) — 1, thereby discarding terms which are 
independent of Xe(* u ) an d Xm(™)i leading to 



x xi(») - - - - 



rl_ ~ b 2 d 2 e- 2bs 



If u 



2 V 6c 



Xm<>) 




u \ 1 u 



be) 2 V 6c 
XcOkmO) j , (B7) 

Xc(™)Xm(m)| , 
(B8) 



where we have set d' — d, Xc(* u ) = Xc(iu), and Xm(* u ) = 
Xm(iu). Substitution of Eqs. (JSTJl and fE5f) into Eq. JTSJl 
leads to Eq. (J7SJ. In order to see that this leading cor- 
rection corresponds to the process of radiation being re- 
flected at the back (left) plate while aquiring finite phase 
shifts upon transmission trough the front (right) plate, 
we note that up to linear order in x' c (^ u )i Xm(* u )j an d 
d' , the terms in square brackets in Eqs. (|B5|I and l|B6|) 
are equal to the phase factor e~ 2h " d , as can be easily 
verified by recalling Eq. I|13fl : 

1 n' 2 (iu)b 2 + b* d , 
/j,'(iu)b 

- 1 - 2b \ 1 + \(j£) [Xe(iu)+Xm(iu)]\d' 



1 - 26 M d' 



-2bLd' 



(B9) 



similar for Eq. JB6I 
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